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Q ' ABSTRACT 

o : 

^SJ ' Context. Recently, the H.E.S.S. Collaboration discovered a very high energy gamma-ray point source close to the Galactic plane. 

They offered three possible explanations for the nature of the source: a gamma-ray binary, a pulsar wind nebula, or a BL Lacertae 

I I ' object. They concluded that the observations favoured an extreme BL Lacertae object interpretation. 

, Aims. We investigated the nature of the radio source reported as the counterpart of the very high energy gamma-ray source. 

a Methods. We performed high-resolution radio interferometric observations with the European Very Long Baseline Interferometry 

' Network (EVN) at a frequency of 1.6 GHz on 2011 May 18. We also reanalysed archival 1.4-GHz radio continuum and HI spectral 

r-j I line data taken with the Very Large Array (VLA). 

Results. The accurate position of the radio source, as observed with EVN, is ~ 4" off from the one obtained in the NRAO VLA 

1 ■ Sky Survey (NVSS). The new position is in excellent agreement with that of the proposed X-ray counterpart of the TeV source. The 
O ' brightness temperature of the compact radio source is ~8 x 10' K. From HI absorption data, a distance of about 1 1.5 ± 1.5 kpc can be 
^ ' inferred for this source. The large-scale HI data unveiled the presence of a shell-like feature with the radio/X-ray/TeV point source in 
C/3 , its interior We interpret this shell as the last vestige of a very old supernova that exploded in a tenuous environment created by the 

■ stellar wind of its massive stellar precursor 

Conclusions. The estimated brightness temperature of the radio point source counterpart of HESS J1943+213 is well below the value 

' expected from the Doppler-boosted radio emission of a BL Lacertae object. This fact and the discovery of traces of a distant supernova 

, explosion around the location of the TeV source lead us to conclude that the most likely origin of the high-energy emission is a remote 

JiT . pulsar wind nebula. If this scenario is true, then the HI shell around HESS J1943-I-213 may represent a population of hitherto missing 

' Galactic supernova remnants. 

' Key words. Techniques: interferometric - Radio continuum: general - Radio lines: ISM - ISM: supernova remnants - Gamma rays: 

lO I general - X-rays: individual: IGR J19443+21 17 

d ■ 

^ ' 1. Introduction ti-a. They concluded that IGR J19443-H2117 is associated with 

' the Chandra source CXOU J194356.2-H21 1823 and the prob- 

I The High Energy Stereoscopic System (H.E.S.S.) consists of ability of spurious association is only 0.39%. According to the 

• r-j four imaging atmospheric Cheren kov telescopes si tuated in the precise Chandra measurement, the X-ray source is located 23" 

^ Khomas Highland of Namibia (lAharonian et al.1 [2006). The away from HE SS J1943-H213 but s till wi thin the H.E.S.S. error 

^ H.E.S.S. Collaboration has been surveying the Galactic plane for circle, leading lAbramowski et alj (1201 lb to identify the X-ray 

new very high energy (VHE >100 GeV) gamma-ray sources, source with the H.E.S.S. source. All the proposed counterparts 

Recently Abramowski etalj (1201 ih reported the discovery of of HESS J1943+213 discussed by Abramowski et al. (2011) ai-e 

an unresolved VHE gamma-ray source close to the Galactic located within the 68 % best-fit source positi on confidence level 

plane, HESS J1943-H213. The conducted search for multi- conto ur of HESS J1943-t-213 (see Fig. 6. in lAbramowski et all 

wavelength counterparts revealed the presence of a hard X- |201 ih . 
ray source, observed by INTEGRAL, IGR J19443-H21 17, in 
the vicinity of the H.E.S.S. source. Landi et al. (2009) ana- 
lyzed the data of the X-ray telescope onboard the Swift satel- In a search for a radio counterpart, a possible source is found 

lite to search for soft X-ray counterparts of three INTEGRAL in the U.S. National Radio Astronomy Observatory (NRAO) 

sources, among them IGR J19443-I-21 17. They concluded that Very Large An^ay (VLA ) Sky Survey (NVS S). The source, 

they found firm soft X-ray localization of the source (SWIFT NVSS J194356-H21 1826 (ICondon etalJll998l) is located 24'.'7 

J1943. 5-1-2120). Later, Tomsick et al. (2009) conducted Chandra away from the H.E.S.S. position (still within the H.E.S.S. enw 

X-ray observations of several INTEGRAL sources including IGR circle) and 3'.'5 away from the Chandra position. The Chandra 

J19443-i-21 17 to localize and measure their soft X-ray spec- source is outside the NVSS error circle. 
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lAbramowski et al.1 (|201 ih discuss possible origins for 
HESS J 1943 +2 13, proposing that it can be either a gamma-ray 
binary, a pulsar wind nebula (PWN), or an active galactic nu- 
cleus (AGN). They conclude that observational facts favour the 
AGN interpretation, and suggest that HESS J 1943 -1-2 13 is an ex- 
treme, high-frequency peaked BL Lacertae object. 

To better understand the nature of the VHE source, we con- 
ducted exploratory Very Long Baseline Interferometry (VLBI) 
continuum observations with the European VLBI Network 
(EVN) at 1.6 GHz frequency of NVSS J 1 94356 -h21 1826 (here- 
after J1943-i-2218). These observations permit to spatially re- 
solve the radio source and to obtain its positional information 
with the precision of a few milli-arcseconds (mas). The high an- 
gular resolution study was complemented with the analysis of 
1 .4-GHz VLA radio continuum and neutral hydrogen emission 
and absorption data in an extended field around the VHE source. 
We describe the observations and data reduction in Sect. [31 and 
discuss our findings in Sect.[3]of this paper Our conclusions are 
summarised in Sect. H] 



2. Observations and data reduction 

The exploratory EVN observation of J1943-H21 18 took place on 
201 1 May 18. At a recording rate of up to 1024 Mbit s"', seven 
antennas participated in this e-VLBI experiment: Effelsberg 
(Germany), Jodrell Bank Lovell Telescope (UK), Medicina 
(Italy), Onsala (Sweden), Toruh (Poland), Hartebeesthoek 
(South Africa), and the phased array of the Westerbork Synthesis 
Radio Telescope (WSRT, The Netherlands). In an e-VLBI exper- 
iment ((Szomoru 2008), the signals received at the remote radio 
telescopes are streamed over optical fibre networks directly to 
the central data processor for real-time correlation. The correla- 
tion took place at the EVN data processor in the Joint Institute 
for VLBI in Europe (JIVE), Dwingeloo, The Netherlands, with 
2 s integration time. The observations lasted for 2 hours. Eight 
intermediate frequency channels (IFs) were used in both right 
and left circular polarisations. The total bandwidth was 128 MHz 
per polarisation. 

The target source was observed in phase-reference mode 
to obtain precise relative positional information. It was cru- 
cial for strengthening the identification of the source, since 
there was a significant difference between the NVSS peak po- 
sition and the Chandra position. The phase-reference calibrator 
source J1946+2300 is separated from the NVSS radio source 
by 1?77 in the sky. Its coordinates in the current 2nd reali- 
sation of the International Celestial Reference Frame (ICRF2) 
are right ascension a ^ 1 9^46'°ff.25 14 0484 and declination 
6 = H-23°0'4'.'4144890 (iFev et al .1120091) . The target-reference 
cycles of ~5min allowed us to spend ~3.5min on the target 
source in each cycle, thus providing almost 1.3 h total integra- 
tion time on J 1943-1-21 18. The source position turned out to be 
offset from the phase centre position (taken from the NVSS cat- 
alogue) by ~4", but still within the undistorted field of view of 
the EVN. 

The NRAO Astronomical Image Processing System (AIPS) 
was used f or the data calibration in the standard way. We refer to 
iFrev et aP (l2008l) for the details of the data reduction and imag- 
ing. The calibrated data were exported to the Caltech Difmap 
package for imaging. Phase self-calibration was only performed 
at the most sensitive antennas (Effelsberg, Jodrell Bank, WSRT). 
No amplitude self-calibration was applied. Finally, the longest 
baselines (from European telescopes to Hartebeesthoek) were 
excluded from the imaging because the signal was barely above 




Right Ascension (mas) 



Fig.l. 1.6-GHz EVN image of J 1943 H-21 18. The lowest con- 
tours are drawn at ±0.6 mJy/beam. The positive contour levels 
increase by a factor of 2. The peak brightness is 25.3 mJy/beam. 
The Gaussian restoring beam shown in the lower-left corner is 
43.9 mas X 28.5 mas (full width at half maximum, FWHM) with 
major axis position angle 40?3. 



the noise level due to the resolved nature of the source. The re- 
sulting image of J1943-I-21 18 is displayed in Fig.[T] 

We also analysed archival VLA continuum data (project 
AH196, observing date 1985 September 30). These observations 
were taken from the region around our target source at 1 .4 GHz, 
in the C configuration of the array, thus improving the angular 
resolution by a factor of three compared to the NVSS survey, 
which was performed with the most compact D configuration of 
the VLA. 

The archival experiment AH 196 included five different 
pointings covering a region of about 2?4xl° around our source 
of interest. The calibration was performed in AIPS, using 3C 48 
as the primary flux density calibrator We used Difmap to obtain 
the image shown in Fig.|2l 

We also investigated the neutral hydrogen (HI) radio emis- 
sion in a large field around the VHE source. To carry out this 
study, we mad e use of the dat a from the VLA Galactic Plane 
Survey (VGPS. IStiTetanilOOel) . 



3. Results and discussion 

3.1. Radio continuum: no blazar-like compact emission 

The phase-referenced exploratory EVN observation pro- 
vided accurate equatorial coordinates for J1943-I-21 18: a = 
19H3'"56!2372 + 0!0001 and 6 = 21°18'23'.'402 + 0'.'002. This 
position agrees well, within uncertainties, with the coordinates 
of CXOUJ194356.2H-211823, the X-ray sou rce proposed to be 
the co unterpart of the H.E.S.S. point source (lAbramowski et all 
Therefore we can confirm that Jl 943 -1-21 18 is indeed 
the radio counterpait of CXOU J194356.2H-211823 and conse- 
quently very likely to be associated with the VHE emission de- 
tected by H.E.S.S. It has to be noted that the EVN detection is 
3'.'75 off from the position given in the NVSS catalogue. 
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Fig. 2. 1.4-GHz VLA C-configuration image of J 1943 +21 18. 
The peak brightness is 60mJy/beam. The lowest contours are 
drawn at + 1 .2 mJy/beam and the positive contour levels increase 
by a factor of 2. The Gaussian restoring beam shown in the 
lower-right corner is 17'.'8 x 15'.' 1 with major axis position angle 
60? 1. Label number 1 indicates the TeV source position and its 
90 % error measured by H.E.S.S. The NVSS source position and 
the position of the X-ray counterpart is labelled as number 2 and 
number 3, respectively. The sizes of the symbols represent the 
eiTors of the corresponding position measurement. The Chandra 
position coincides with the position of the radio source derived 
from our phase-referencing EVN observation. (The accuracy of 
the EVN position is superior to that of the X-ray, therefore it 
cannot be distinctively displayed in this figure.) 



We used the Difmap package to fit a circular Gaussian 
brightness distribution model component to the VLBI visibil- 
ity data. The feature in our EVN image can be well described 
with a component of 31 mJy flux density and 15.8 mas angular 
size (full width at half maximum, FWHM). These values imply 
a brightness temperature of Tb = 7.9xlO^K. Since Jl 943 +21 18 
is very close to the Galactic plane (at about -1°.3 Galactic lati- 
tude), angular broadening caused by the intervening ionised in- 
terstellar matter can distort the image of a dista nt compact extra- 
alactic radio source. According to the model of lCordes & Laziol 
200l . the maximal amount of angular broadening of a point 
source in this direction at 1 .6 GHz (at the frequency of our EVN 
observation) is expected to be 3.34 mas. However, even if we 
take this effect into account, the "de-broadened" source size re- 
mains quite large. The brightness temperature calculated using 
this "de-broadened" source size is still substantially lower than 
the intrinsic equipartit ion value estimated for relativistic com- 
pact jets (~ 5 X 10'° K, lReadheadl994) . Thus relativistic beam- 
ing most probably does not play a role in the appearance of this 
source. 

The flux density recovered in our high-resolution EVN ob- 
servation is only one-third of the value reported by NVSS at 
1.4 GHz, (102.6±3.6)mJy. To investigate the discrepancy be- 
tween the flux density values, we analysed the WSRT phased 
aiTay data taken during our EVN experiment. The offset of the 
source position from the phase centre is a significant fraction 
of the synthesised beam of the phased array of the WSRT (7'.'8 
at 1 .6 GHz) therefore the derived source characteristics should 




Fig. 3. Radio continuum emission at 1 .4 GHz as taken from the 
VGPS (IStil et al. 2006). The white cross indicates the location 
of J 1943+21 18. Note that this map is presented in galactic coor- 
dinates. 



be used with caution. However, the obtained flux density value, 
~95 mJy, still agrees well with the value reported in the NVSS 
for J 1943+21 18. Thus, we can conclude that the high-resolution 
EVN observation resolved out significant portion of the large- 
scale structure of J 1943+2 118. This was confirmed with the 
archival VLA C-array data (see Fig. [J]). According to those, the 
flux density of the source is (91±5)mJy, which agrees with the 
flux density values derived from the two other (the NVSS and 
the WSRT-only) lower-resolution data sets. 

Additionally, the VLA C-array data revealed an elongated 
structure of Jl 943 +21 18 in the SE-NW direction. This shape can 
naturally explain why the NVSS position is ~ 4" off. In Fig.|2] 
label 2 indicates the NVSS position, while label 3 indicates the 
coinciding Chandra X-ray and EVN radio positions. The shift 
between points 2 and 3 is exactly in the direction of the source 
extension. 

Finally, to gain insight into the field where HESS J1943+213 
is located, and to search for possibly associated extended emis- 
sion, we inspected the VGPS image of the radio continuum 
emission at 1.4 GHz in a 5°x5° area around the source posi- 
tion (Fig.|3]l. Essentially, it shows an almost empty field around 
J1943+21 18 (whose location is indicated by a white cross in Fig. 
|3]l, with no trace of diffuse emission at the sensitivity of this 
survey. The reprocessed VLA C-aiTay archival data confirm this 
picture, down to an rms noise level of 0. 15 mJy/beam. 

3.2. HI data: the distance and the associated features 

With an angular resolution of 1' and a sensitivity of 1 1 mJy/beam 
(2K for a 1' beam), the VGPS (I Stil et al.ll2006h is an adequate 
database to search for possible footprints of the event that re- 
sulted in the very high-energy emission detected by H.E.S.S. We 
have carried out an HI absorption study to constrain the distance 
to J1943+21 18. Fig.|4]displays the HI profile traced towards the 
point source (top panel) and the absorption spectrum obtained 
from the subtraction of this profile from the average of four sur- 
rounding offset points (lower panel). From these spectra we can 
conclude that there is an absorption feature at ya;- 16 km s"' , fol- 
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Fig. 4. Top panel: HI spectrum towards J1943+2118. Bottom 
panel: absorption HI spectrum obtained from the subtraction of 
the above profile from the average of four offset profiles. 



lowed by an emission feature peaking near c;^-50kms ~ . By 
applying a Galactic circular rotation model (Fich et al. 1989|). 
we set tentative limits for the distance to J1943+21 18, conclud- 
ing that it is beyond lOkpc (corresponding to the radial velocity 
-16km s ') but closer than 13kpc (corresponding to the radial 
velocity -50km s"'). In what follows we adopt as a compro- 
mise an intermediate distance of 11.5kpc for J1943-I-21 18. We 
note that the emission feature seen at -50 km s ~' lies most prob- 
ably i n the warping region of the outer galaxy (IVoskes & BurtonI 
120061) further than Jl 943 -1-21 18 and does not seem to be related 
to it. 

By inspecting the HI emission in the large 5°x5° field, across 
the whole observed velocity range between -114kms"' and 
+ 166 km s ' , we find that in all channel maps between radial ve- 
locity i;a;-H50kms"' and (;a;-i-57kms"', a striking, almost com- 
plete shell-like feature surrounds the location of Jl 943 -1-21 18 
(Fig.H). 




Fig. 5. HI emission distribution from VGPS, within the velocity 
range -H50kms"' to -i-57kms"'. The white cross shows the lo- 
cation of J1943-)-21 18 and the white dashed circle marks up the 
location of the discovered HI shell. Note that this map is pre- 
sented in galactic coordinates. 




Fig. 6. A detailed view of the HI shell shown in Fig.|5] The cross 
shows the location of J1943-I-21 18 and the circle is the 2f8 con- 
fidence size of HESS J1943+213. 

This HI shell, displayed in detail in Fig. |6] appears at 
a "forbidden" velocity for the expected rotational properties 
of the Galactic gaseous disk, since in this direction of the 
Galaxy (Galactic longitude l~57°) the maximum positive veloc- 
ity (corresponding to the tangent point) is predicted to be about 
+35kms This fact is not surprising. Recently iKang et aH 
discussed the presence of many faint, wing-like HI fea- 
tures at velocities beyond the boundaries allowed by Galactic 
rotation (what they call Forbidden- Velocity Wings, FVWs), and 
conclude that these features are the results of violent events, 
probably missing old Galactic supernova remnants (SNRs) 
which are invisible in other spectral regimes. 

With this idea in mind, we investigated the characteris- 
tic physical parameters of the newly detected HI shell around 
J1943-(-21 18. In Table[T]we summarise the observed and derived 
parameters of the HI shell (for the adopted distance of 1 1 .5 kpc). 
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Table 1. Observed and derived parameters of the HI shell. 



Geometrical centre: 




Equatorial coordinates 


a = 19H3'"23%(5 = +21° 10' 45" 


Galactic coordinates 


! = 57°.6,b = -1°.25 


Angular diameter 


r 


Linear diameter 


200 pc 


Adopted systemic velocity 


-16kms"' 


Adopted central velocity 


+54kms"' 


Expansion velocity 


70kms"' 


HI mass 


6.4 X 10"* Mo 


Ambient density 


0.6cm-3 


Kinetic energy 


1.6 X lO^'erg 


Initial energy 


1.6 X 10^2 erg 


Dynamic age 


4 X 10' years 



The central velocity was derived on the basis of the velocity of 
the central channel of the nine where the HI shell is clearly de- 
tected between +50kms"' and +57kms"'. The expansion ve- 
locity is estimated by adding the ~54km s~' at which the shell is 
seen, to the absolute value of the assumed systemic velocity of 
approximately -16kms"', as suggested by the absorption study. 
We also searched for the expanding "caps" of the shell at the ap- 
propriate velocities, i. e. (-16 - 70)kms"'= -86kms"'. Even 
though quite faint, we found traces of the HI shell at a very close 
velocity value, -81 kms"'. The HI mass was calculated by inte- 
grating all contributions within the shell across all the channels 
where the shell is present. The ambient density was derived by 
assuming that all the swept-up mass that forms the shell was ini- 
tially uniformly distributed within a sphere of 200 pc diameter. 
The initial energy was calculated from the relation (Ciofh et al. 
[1988 ): £sN = 6.8 x lO"*^ n^"^ R^^^ v^^^^ ^ i^' ergs, where ^, 
the metallicity index, is a dopted 0.2, corresponding to t/^ilOkpc 
dMaciel & Ouirezalll999h . Finally, the dynamic age is calculated 
from the relation f ~ 0.37?s/uexp- 

All the observed and calculated paramete rs for this shell ar e 
completely analogous to those obtained by iKoo et all (l2006h . 
who conducted an HI study towards one of the FVWs, the struc- 
ture FVW 190.2-1-1.1. They concluded that such parameters are 
only consistent with the HI shell being the last remnant of a su- 
pernova (SN) explosion that occurred in the outermost fringes 
of the Galaxy some 3x10^ years ago, an SNR that is not seen 
in any other wavebands, and that represents the oldest type of 
SNRs, essentially invisible except via its HI line emission. 

In our case, the observed HI shell might be the last vestige of 
an old, very energetic SN that exploded on the other side of the 
Galaxy, about 400 000 years ago, and is currently expanding in a 
tenuous interstellar gas. These facts explain the absence of radio 
continuum emission associated with this remnant. In this sce- 
nario, the sources detected in TeV, radio, and X-rays would be 
the emission from the PWN. For comparison, if the paradigmatic 
PWNe Crab and 3C 58 were located at a distance of Il.Skpc, 
as the approximate distance of J 1943 -1-2 11 8, their angular size 
in radio wavelengths would be of the order of 1'2 x 0.'8 for 
Crab and ~0.'8 for 3C 58, in excellent concordance with the 
size of the elliptical structure of J1943-H2118 as seen from the 
VLA observation. The size of the PWNe in X-rays is generally 
smaller than in radio due to the smaller synchrotron lifetimes 
of the higher-energy electrons (Slane et al. 2000), thus explain- 
ing why Chandra detected it as a point -like source. Besides , the 
radio spectral index -0.32 reported bv lVoIImer et al.l (l2OI0l) for 



the NVSS point source coincident with J 1943 -1-2 11 8 is compat- 
ible with the standard spectral indices of PWNe (between and 
-0.3, [Gaensler & Slane 2006). 

However, one p roblem arises in this scenario. According to 
iMattana et al.l (12009) . the gamma-ray to X-ray flux ratio is pro- 
portional to the age of the PWN. The relatively low value of the 
flux ratio (0.04) of J 1943 -1-21 18 implies that the nebula is only 
~10^ years old, more than hundred times younger than the age 
implied by the dynamics of the HI shell. This contradiction can 
be resolved if the HI shell was produced by the stellar wind of the 
massive supernova progenitor. As the SNR expanded in a previ- 
ously evacuated cavity, the ambient density to create the SNR 
shell was very low and it is natural that the SNR remained un- 
detected. Such HI bubbles are common features around massive 
early-type stars which eventually end their lives e xploding as 
SNe (e.g. Dubneretal. 1990; CacDaetal. 1996; Va squez et aP 
I2OO5I: iGiacani etaIJ I2011). In thiV scenario, the estimated dy- 
namic age has to be attribut ed to the wind-blown shell. A similar 
scenario was suggested by iKoo & HeilesI (Il995h to explain the 
detection of a high-velocity HI shell associated with the SNR 
W44. 

4. Conclusions 

We reported on the study of the radio counterpart of the newly 
discovered TeV source, HESS J1943-I-213. We analysed new 
and archival radio continuum data and archival HI spectral- 
line observati ons of the pr oposed radio counterpart, J 1943 +21 18 
( Abramowski etalJl20Ilh . 

Our high-resolution exploratory EVN observation provided 
the most accurate position of the source to date, thus strength- 
ening the identification. From the EVN data, we could also cal- 
culate the brightness temperature of the source. The resulting 
low value indicates that relativistic beaming most probably does 
not play a role in the appearance of t his source. According t o 
the unified model of radio-loud AGNs (IUrrv&PadovanillI995h . 
however, BL Lacertae objects are viewed at angles very close 
to the jet direction, therefore their radio emission is relativis- 
tically beamed, with apparent brightness temperatures well in 
excess of the equipartition value. Therefore, the measured low 
brightness temperature of J 1943 +22 18 is not compatible with 
the BL Lacertae object n ature, thus weakening the suggestion of 
lAbramowski et alJ (12011 ). 

The re-analysed archival VLA C-array data, which recover 
emission at spatial scales that have been filtered out in the EVN 
observation, revealed an elongated source structure. Further ra- 
dio observation at intermediate resolution is needed to better un- 
derstand the radio structure of J 1943+22 18. 

On the other hand, the HI data revealed a shell-like feature 
around JI943+2II8. Even though we do not have yet an un- 
equivocal evidence confirming that the two objects are associ- 
ated, all the observed and derived parameters of this large-scale 
feature and the properties of the radio continuum source point to 
the discovery of a distant SNR whose PWN could be responsible 
for the TeV emission detected by H.E.S.S. If the HI shell is ex- 
plained as a result of a SN explosion, then its dynamics indicate 
that the SN explosion took place ~ 400 000 years ago. The X- 
ray to gamma-ray flux ratio on the other hand suggests a much 
younger (~ 10^ yrs) PWN. If the latter age is the correct, then the 
observed HI shell was not produced by the SN explosion but by 
the action of the stellar wind of the massive SN progenitor(s). 
Later the SN explosion took place within the rarefied wind cav- 
ity and the nebula created by its pulsar gave rise to the observed 
gamma-ray and X-ray emission. An ultimate evidence for the 
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PWN scenario would be the detection of pulsed radio emission 
from J1943+21 18. It is worth mentioning that the SNR W44, a 
high energy gamma-ray emitter with an associated high-velocity 
HI shell, also harbours a pulsm" (PSR B1853+01) which is pow- 
ering a PWN jFrail et alil9"96l) . 

The importance of the reported discovery resides not only in 
the fact that it naturally explains the nature of the TeV emission, 
but also, as remarked by Koo et al. (2006), it helps to solve the 
long-standing problem of the "missing" Galactic SNRs, where 
the detected SNRs barely make up 1 % of the expected number. 
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